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Abstract 

We consider the well-posedness of a model for a flow-structure interaction. This model 
describes the dynamics of an elastic flexible plate with clamped boundary conditions im- 
mersed in a supersonic flow. A perturbed wave equation describes the flow potential. The 
plate's out-of-plane displacement can be modeled by various nonlinear plate equations (in- 
cluding von Karman and Berger). We show that the linearized model is well-posed on the 
state space (as given by finite energy considerations) and generates a strongly continuous 
semigroup. We make use of these results to conclude global-in-time well-posedness for the 
fully nonlinear model. 

The proof of generation has two novel features, namely: (1) we introduce a new flow po- 
tential velocity-type variable which makes it possible to cover both subsonic and supersonic 
cases, and to split the dynamics generating operator into a skew-adjoint component and a 
perturbation acting outside of the state space. Performing semigroup analysis also requires 
a nontrivial approximation of the domain of the generator. And (2) we make critical use 
of hidden regularity for the flow component of the model (in the abstract setup for the 
semigroup problem) which allows us run a fixed point argument and eventually conclude 
well-posedness. This well-posedness result for supersonic flows (in the absence of rotational 
inertia) has been hereto open. The use of semigroup methods to obtain well-posedness 
opens this model to long-time behavior considerations. 

Key terms: flow-structure interaction, nonlinear plate, supersonic and subsonic flows, 
nonlinear semigroups, well-posedness, dynamical systems. 

MSC 2010: 35L20, 74F10, 35Q74, 76J20 



Introduction 



1.1 Physical Motivation 

The interaction of a thin, flexible structure with a surrounding flow of gas is one of the principal 
problems in aeroelasticity. Models of this type arise in many engineering applications such as 
studies of bridges and buildings in response to wind, snoring and sleep apnea in the human palate, 
and in the stability and control of wings and aircraft structures [TJ [5J HU [221 133] • In general, for 
an abstract setup, we aim to model the oscillations of a thin flexible structure interacting with 
an inviscid potential flow in which it is immersed. These models accommodate certain physical 
parameters, but one of the key parameters is the flow velocity of the unperturbed flow of gas. 

Specifically, we deal with a common flow-structure PDE model which describes the interac- 
tive dynamics between a (nonlinear) plate and the surrounding potential flow (see, e.g., [71 122]). 
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This model is one of the standard models in the applied mathematics literature for the modeling 
of flow-structure interactions (see, e.g., [71(22] and also [201 [H| and the references therein). 

The main goal of this paper is to present Hadamard well-poscdness results for the model in 
the presence of supersonic flow velocities. While subsonic flows have received recent attention 
which has resulted in a rather complete mathematical theory of well-posedness [TUl HSJ [UJ 03] 
and spectral behavior for reduced (linear) models (3] [39] , this is not the case for the supersonic 
flow velocitie^. The mathematical difficulty in going from a subsonic to supersonic regimes is 
apparent when one inspects the formal energy balance. There is an apparent loss of ellipticity 
affecting the static problem. This, in turn, leads to the appearance of boundary trace terms that 
can not be handled by known (elliptic) PDE-trace theories. Successful handling of this issue 
yields new methodology which is based on appropriate (microlocal) boundary trace estimates 
and effectively compensates for this loss of ellipticity. The method here presented additionally 
covers (with minimal adjustments) subsonic flows. 

Thus, with respect to the supersonic model, this paper addresses the open question of well- 
poscdness of finite energy solutions, which is the most fundamental for future studies. Well- 
poscdness results are necessary mathematically in order to begin long-time behavior and control 
studies of the model, which belong to the most interesting and pertinent mathematical studies 
in application for PDE models. Having shown well-posedness allows us to move into stability 
studies in the presence of control mechanisms [35J ISO] • 

1.2 Notation 

For the remainder of the text we write x for (x, y,z) G R+ or (x, y) G f2 C ffi.?^ y \y, as dictated by 
context. Norms || • || are taken to be L2(D) for the domain dictated by context. Inner products 
in i2(R+) are written (•, •), while inner products in L2(M 2 = <9R+) are written <■,■>. Also, 
H S (D) will denote the Sobolev space of order s, defined on a domain D, and Hq(D) denotes 
the closure of C£°(D) in the H S (D) norm which we denote by || • \\h s (D) or II ' IIs.d- We make 
use of the standard notation for the trace of functions defined on , i.e. for <f> £ H 1 (M.^_), 
j[4>] = a _ is the trace of cj> on the plane {x : z = 0}. 

1.3 PDE Description of the Model 

The model in consideration describes the interaction between a nonlinear plate with a field or 
flow of gas above it. To describe the behavior of the gas we make use of the theory of potential 
flows (see, e.g., [7j [20] and the references therein) which produce a perturbed wave equation 
for the velocity potential of the flow. The oscillatory behavior of the plate is governed by the 
second order (in time) Kirchoff plate equation with a general nonlincarity. We will consider 
certain 'physical' nonlinearities which are used in the modeling of the large oscillations of thin, 
flexible plates - so-called large deflection theory. 

The environment we consider is = {(x,y,z) : z > 0}. The plate is modeled by a 
bounded domain fl C = {(x, y, z) : z = 0} with smooth boundary dfl = T. The 

plate is embedded in a 'large' rigid body (producing the so-called clamped boundary conditions) 
immersed in an inviscid flow (over body) with velocity U =/= 1 in the negative a;-directior|^|. This 
situation corresponds to the dynamics of a panel element of an aircraft flying with the speed U, 
see, e.g., [2"T] . 

J We exclude the models which are based on the so-called "piston" theory, see [7] Chapter 4], 1221 Part I], and 
also Remark 6.2.2 in 1161 for a recent discussion. 

2 Here we normalize U = 1 to be Mach 1, i.e. < U < 1 is subsonic and U > 1 is supersonic. 
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The scalar function u : x R + — > K represents the vertical displacement of the plate in the 
z-direction at the point (x; y) at the moment t. We take the nonlinear Kirchoff type plate with 
clamped boundary condition^: 

' u t t + A 2 u + f(u) = p(x,t) inOx(0,T), 

u(0) = u ; u t (0)=ui, (1.1) 
u = d„u = onafix(0,T). 

The aerodynamical pressure p(x, t) represents the coupling with the flow and will be given below. 

In this paper we consider a general situation that covers typical nonlinear (cubic-type) force 
terms f(u) resulting from aeroelasticity modeling [7J HOI HD HI] • These include: 

Assumption 1.1. 1. Kirchhoff model: u n- f(u) is the Nemytski operator with a function 
/ £ Lipi oc (R) which fulfills the condition 

f(s) 

liminf^- > -Ai, (1.2) 

|s|->oo S 

where Ai is the first eigenvalue of the biharmonic operator with homogeneous Dirichlet 
boundary conditions. 

2. Von Karman model: f(u) = —[u, v(u) + Fq], where Fo is a given function from H 4 (£l) and 
the von Karman bracket [u, v] is given by 

[«, v] = d 2 x u ■ d 2 y v + d 2 y u ■ 8 2 v - 2 • d 2 xy u ■ d 2 y v, 

and the Airy stress function v(u) solves the following elliptic problem 

A 2 v(u) + [u,u] = in CI, d v v(u) = v(u) = on dil. (1.3) 

Von Karman equations are well known in nonlinear elasticity and constitute a basic model 
describing nonlinear oscillations of a plate accounting for large displacements, see [26] and 
also [HI HI] and references therein. 

3. Berger Model: In this case the feedback force / has the form 



/(«) 



\Vu\ 2 dx-T 



Am, 



where K > and T £ E are parameters, for some details and references see [5] and also 
[131 Chap.4]. 

For the flow component of the model, we make use of linearized potential theory, and we 
know [SI [7] [3T] that the (perturbed) flow potential cf> : — > R must satisfy the perturbed wave 
equation below (note that when U — this is the standard wave equation): 

\d t + Ud x ) 2 = A(b inR^x(0,T), 

0(0) = O ; M0) = 4>i, (1-4) 
d u <f> = d(x,t) onR| M} x(0,T). 



3 While being the most physically relevant boundary conditions for the flow-plate model, clamped boundary 
conditions allow us to avoid certain technical issues in the consideration and streamline our exposition. Other 
possible and physically pertinent plate boundary conditions in this setup include: hinged, hinged dissipation, 
and combinations thereof 1271 . 
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The strong coupling here takes place in the downwash term of the flow potential (the Neumann 
boundary condition) by taking 

d(x,t) = -[(d t + Ud x )u(x)] .1q(x) 

and by taking the aerodynamical pressure of the form 

p{x,t) = {d t + Ud x ) 1 [<t>] (1.5) 

in above. This gives the fully coupled model: 

'u tt + A 2 u + f(u) = (d t + Ud x )~/[(f>\ infi x (0,T), 
u(0) = u o ; ut(0)=ui, 

u = d v u = ondQx(0,T), 
(9 t + Ud x ) 2 4> = A0 inK'|x(0,T), 
#O)=0oi MO) = <h, 

d„(t> = -[{d t + Ud x )u(x:)] .l n (x) on K {(^ y) } x (0,T). 
1.3.1 Parameters and New Challenges 

We do not include the full rotational inertia term M a = (1 — aA)u tt in the LHS of plate equa- 
tion, i.e., we take a = 0. In some considerations (see [27]) this term is taken to be proportional 
to the cube of the thickness of the plate, however, it is often neglected in large deflection theory. 
From the mathematical point of view, this term is regularizing in that it provides additional 
smoothness for the plate velocity u t , i.e. £2^) — > H 1 ^). This is a key mathematical assump- 
tion in our analysis which separates it from previous supersonic considerations and increases 
the difficulty of the analysis. This will be further elaborated upon in the discussion of previous 
literature below. The case a > presents modeling difficulties in problems with flow coupling 
interface, but is often considered as a preliminary step in the study limiting problems as a \ 
(see OHO], [31], and also [16]) where subsonic regimes were studied. On the other hand, the 
case when rotational terms are not included (a = 0) leads to substantial new mathematical 
difficulties due to the presence of flow trace terms interacting with the plate. 

The second key parameter is the unperturbed flow velocity U. Here we take U 7^ 1 arbitrary. 
However, the supersonic case (U > 1) is the most interesting case from the point of view of ap- 
plication and engineering. Results in this case can be more challenging, due to the loss of strong 
cllipticity of the spatial flow operator in (|1.4p . For the subsonic case (0 < U < 1) there are other 
methods available, see, e.g., [8] [16j [17l [43] . However, for the non-rotational case a = in the 
supersonic regime U > 1, the problem of well-posedness of finite energy solutions is challenging 
and has been hereto open. As is later expounded upon, the lack of sufficient differentiability and 
compactness for the plate velocity component u t renders the existing methods (see |16( Sections 
6.5 and 6.6], for instance) inapplicable. 

The aim of this paper is to provide an affirmative answer to the well-poscdncss question in 
the (mathematically) most demanding case with a = and U > 1. In fact, we will show that 
the resulting dynamics generate a nonlinear semigroup associated with mild solutions. Though 
in this treatment we focus on the most challenging case: a = and U > 1, the new methods 
developed apply to the full range U 7^ 1. 

1.4 Energies and State Space 

In the subsonic case < U < 1 energies can be derived by applying standard multipliers ut and 
4>t along with boundary conditions to obtain the energy relations for the plate and the flow. 
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This procedure leads to the energy which is bounded from below in the subsonic case. However, 
it is apparent in the supersonic case that we will obtain an unbounded (from below) energy of 
the flow. Hence, we instead make use of the flow acceleration multiplier (dt + Ud x )4> = ip. Our 
so-called change of variable is then <j> t — y (<j) t + U<f> x ) = ip. Thus for the flow dynamics, instead 
of (6 \4>t) we introduce the phase variables (</>; -0). 

We then have a new description of our coupled system as follows: 

'(d t + ud x )4> = ip 
(d t + Ud x )i) = A<t> 
< d u cj> = -[(d t + Ud x )u(x.)] -lo(x) 
utt + A 2 u + f(u) = 
u = d u ii = 

This leads to the following (formal) energies, arrived at via Green's Theorem: 

^w4(lhll 2 + HA u || 2 )+n( u ), (i-8) 

E f i(t)=l(\\n 2 +\\vct>\\ 2 ), 

£{t) =E vl (i) + E fl (t), 

where Tl(u) is a potential of the nonlinear force f(u), i.e. we assume that f(u) is a Frechet 
derivative of H(u), /(it) = W(u). Hypotheses concerning H(u) are motivated by the examples 
described in Assumption 11.11 and will be given later (see the statement of Theorem 13.101) . 
With these energies, we have the formal energy relatiorQ 

£(t) + U f < u x ,7|$ > dt = £(0). (1.9) 
Jo 

This energy relation provides the first motivation for viewing the dynamics (under our change 
of phase variable) as comprised of a generating piece and a perturbation. 

Finite energy constraints manifest themselves in the natural requirements on the functions 
4> and u: 

0(x,t) G C(0,T;i7 1 (R3))nC 1 (0,T;i 2 (]R3)) i (i. 10 ) 
u(x,t) € C(0,T;ff 2 (r!))nC7 1 (0,T; J L 2 (fi)). (1.11) 

In working with well-poscdness considerations (and thus dynamical systems), the above finite 
energy constraints lead to the so-called finite energy space, which we will take as our state space: 

Y = Y f i x Y pl = (H^Rl) x L a (R+)) x (flo(0) x L 2 (fl)). (1.12) 

(2) 

Remark 1.1. The energy Efi(t) defined above coincides with the energy (t) introduced in 
[9l[T0], see also [HI Sect. 6. 6]. As previously indicated, in the subsonic case, the standard flow 
multiplier <j) t is used in the analysis, rather than %j) = <j) t + U<j) x ; this produces differing flow 
energies and an interactive term Ei nt (t), which does not appear here (see, e.g., [51 [TrJl [T71 
More specifically, the flow component of the energy in this case is given by 

E$(t)^±{\\M 2 + \\v<t>\\ 2 - u 2 \\d x <j>\\ 2 ) 

4 For some details in the rotational inertia case we refer to |10| . see also the proof of relation (6.6.4) in [161 
Section 6.6]. 



in R3 x (0,T), 
in R3_ x (0,T), 

onR{ M} x(0,T), (1.7) 
infix (0,T), 
on dn x (0,T). 
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and the interactive energy Ei n t = U < ~f[<f>], d x u >. The total energy defined as a sum of the 
three components £(t) = E^\t) + E p i(t) + E int (t) satisfies £(t) = £(s). We note, that in the 

supersonic case the flow part of the energy E$(t) is no longer nonnegative. This, being the 
source of major mathematical difficulties, necessitates a different approach. In fact, the new 
representation of the energies as in (|1.8p provides good topological measure for the sought after 
solution, however the energy balance is lost in (|1.9[) and, in addition, the boundary term is 
"leaking energy" and involves the traces of L 2 solutions, which are possibly not defined at all. 

In view of the above, our strategy will be based on (i) developing theory for the traces 
of the flow solutions; (ii) counteracting the loss of energy balance relation. The first task 
will be accomplished by exploiting sharp trace regularity in hyperbolic Neumann solutions (see 
[2H EH EH H2] for related results). The second task will benefit critically from the presence of 
the nonlinearity. 

1.5 Definitions of Solutions 

In the discussion below, we will encounter strong (classical), generalized (mild), and weak (vari- 
ational) solutions. In our analysis we will be making use of semigroup theory, hence we will work 
with generalized solutions; these are strong limits of strong solutions. These solutions satisfy 
an integral formulation of (|1.6[) . and are called mild by some authors. In our treatment, we will 
produce a unique generalized solution, and this, in turn, produces a unique weak solution, see, 
e.g., QI1 Section 6.5.5] and g3]. 

We now define strong and generalized solutions: 

Strong Solutions. A pair of functions ((/>(x, y, z;t);u(x,y;t)) satisfying (|1.10|) and (jl.lip is 
said to be a strong solution to (|1.6[) on [0, T] if 

• (font) G i 1 («,&;^ 1 (R+) x H$(fl)) and € L 1 (a, b; L 2 ($t+) x L 2 (fi)) for any 
[a,b] C (0,T). 

• A 2 u(t) -U^[d x (f>(t)] G L 2 (fi) (thusu(t) 6 H 7 / 2 (n)r\H$(n)) and the equation u tt + A 2 u + 
f(u) = p(x,i) holds in H' 1 / 2 ^) for t G (0,T) with p given by (fl~5j) . 

• (U 2 - l)d 2 (f>{t) - (d 2 + d 2 z )4>(t) G L 2 (Rl) with boundary conditions d v <t>(t) G H 1 ^ 2 ) for 
all t G (0,T) and satisfying the relation d v <j> = - [(St + [/3 x )w(x)] ■ l n (x) on R 2 x (0,T). 
Moreover (9 t + Ud x ) 2 (j) = A<f) holds for almost all t G (0, T) and (a;, y, z) G K^_. 

• The initial conditions are satisfied: 0(0) = 0t(O) = 0i, u(0) = uq, u t (0) = U\. 

Remark 1.2. The smoothness properties in the definition above are motivated by the description 
of the generator of the linear problem in the supersonic case U > 1 which is given below, see 
relation (|2.2[) and Lemma 12.11 In the subsonic case regular solutions display more regularity 
(see, e.g., [Ml Sccions 6.4 and 6.5] and [E]). The above analysis also reveals that the degraded 
differentiability of strong solutions is due to the the loss of clipticity in the supersonic regime 
and non-Lopatinski character of the boundary conditions. 

As stated above, generalized solutions are strong limits of strong solutions; these solutions 
correspond to semigroup solutions for an initial datum outside of the domain of the generator. 

Generalized Solutions. A pair of functions (<f>(x, y, z; t); u(x, y; t)) is said to be a generalized 
solution to problem (|1.6[) on the interval [0, T] if (|1.11[) and (|1.10[) arc satisfied and there exists 
a sequence of strong solutions (<j> n (t); u n (t)) with some initial data (0q , ; Uq ; u") such that 

lim max \\\d t (j) - d t <t> n {t)\\ L ^ a ) + \\<f>{t) - to{t)\\m(R 3 ) \ = 
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lim max \\\d t u(t) - d t u n (t)\\ L ^ n) + \\u(t) - u n (t)\\ H 2 {n} \ =0. 

We can show that in the case when the nonlinear term / is locally Lipschitz from Hq(CI) into 
L 2 {£1) the generalized solutions are in fact weak solutions, i.e., they satisfy the corresponding 
variational forms (see Definition 6.4.3 in [T|)J Chapter 6]). This can be verified for strong solutions 
by straightforward integration with the use of regularity exhibited by strong solutions. Using 
the (strong) limit definition of generalized solutions, we can pass in the limit and show that the 
generalized solution satisfies the weak formulation of (|1.6|) . This weak solution is in fact unique 
- for the proof, we defer to the method presented in |16[ Chapter 6]. 



1.6 Description of Past Results 

Flow-structure models have attracted considerable attention in the past mathematical literature, 
see, e.g., [21 EJ HI IE HH1 [HI [HI [HI I22J EH [3H ESJ EZl EH HH S3] and the references therein. 
However, the vast majority of the work done have been devoted to numerical and experimental 
studies, see, e..g., [TJ HJ [7J HOI IHl I2H HI] an d also the survey [33] and the literature cited there. 
Much of the studies has been based on linear one-dimcnsional-spccial geometries plate models 
where the goal was to determine the speed at which flutter occurs, see [1] [6] [7J [21] [24] [34] for 
instance. More recently the study of linear models with a one dimensional structure (beam) 
and Kutta-Jukovsky boundary conditions found renewed interest and have been extensively 
pursued in [2] [39] [40] [41]. This line of work has focused on spectral properties of the system, 
with particular emphasis on identifying aeroelastic eigenmodes corresponding to the associated 
Possio integral equation. 

In contrast, our interest here concerns PDE aspects of the problem, including the funda- 
mental issue of well-poscdncss of finite energy solutions corresponding to nonlinear flow-plate 
interaction in the principal case for the parameters a and U with clamped plate boundary 
conditions. 

In all parameter cases, one is faced with low regularity of boundary traces due to the failure of 
Lopatinski conditions (unlike the Dirichlet case [38], where there is no loss of regularity to wave 
solutions in their boundary traces). In fact, the first contribution to the mathematical analysis 
of the problem is 0[TU] (see also [H] Section 6.6]), where the case a > is fully treated. The 
method employed in [5J [TOJ [TB] relies on the following main ingredients: (1) sharp microlocal 
estimates for the solution to the wave equation driven by iJ 1 / 2 (f2) Neumann boundary data 
given by u t + Uu x . This gives </> t |n € £2(0, T; 1 / 2 (S1) [33] (in fact more regularity is presently 
known: iJ~ 1 / 3 (il) S3); and (2) the regularizing effects on the velocity u t (i.e. u t € 7J 1 (fi)) 
rather than just L2(Cf)), when a > 0. The above ingredients, along with an explicit solver 
for the 3-dimensional wave equation and a Galcrkin approximation for the structure allows one 
to construct a fixed point for the appropriate solution map. The method works well in both 
cases < U < 1 and U > 1. Similar ideas were used more recently [36] [37] in the case when 
thermoelastic effects are added to the model; in this case the dynamics also exhibit _ff 1 (f2) 
regularity of the velocity in both the rotational and non-rotational cases due to the analytic 
regularizing effects induced by thermoelasticity |30j . However, when a = 0, and thermoelastic 
smoothing is not accounted for, there is no additional regularity of Ut beyond Z>2(Q). In that 
case the corresponding estimates become singular. This destroys the applicability of previous 
methods. In summary, much of the work on this problem to date has assumed the presence 
of additional regularizing terms in the plate equation, or depends critically on the condition 
U < 1. A recent book ([TB] Chapter 6]) provides an account of relevant results, including more 
recent applications of the compactness method in the case a — and < U < 1. Existence of 
a nonlinear semigroup capturing finite energy solutions has been shown in [43] . see also [17) . 
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In this treatment we take a = 0, corresponding to the more difficult non-rotational model, 
and we approach the problem with < U =/= 1 from the semigroup point of view - without any 
reliance on explicit solvers for the flow equation or Galerkin constructions. The advantage of this 
approach, in addition to solving the fundamental wcll-posedncss question, is the potential for 
an array of important generalizations, including more general flow equations and more general 
nonlinearities appearing in the structure. Moreover, it may be possible to view the supersonic 
global solution we arrive at for a = as the uniform limit of solutions as a j (as in the 
subsonic case [3"T]). 

The main mathematical difficulty of the problem under consideration is the presence of the 
boundary terms: ((f>t + U4> x ) |n acting as the aerodynamical pressure in the model. When U — 0, 
the corresponding boundary terms exhibit monotone behavior with respect to the energy inner 
product (see [11], [16j Section 6.2] and [28]) which is topologically equivalent to the topology 
of the space Y given by (|1.12[) . The latter enables the use of monotone operator theory ( [TTJ , 
[TBI Section 6.2] and [28]). However, when U > this is no longer true with respect to the 
topology induced by the energy spaces. The lack of the natural dissipativity for both interface 
traces, as well as the nonlinear term in the plate equation, make the task of proving well- 
poscdness challenging. In the subsonic case, semigroup methods were applied to the problem 
by implementing certain bounded adjustments to the inner product structure of the state space 
which then produced shifted dissipativity |43) . This type of consideration is not possible here, 
owing to the degeneracy of the standard energy functional when U > 1 . 

In contrast to these works, our method and results do not depend on any smoothing mech- 
anism (as we take a = 0) and are applicable for all U =/= 1. The key ingredients rely on the 
development of a suitable trace theory for the velocity of the flow and implementing the cor- 
responding estimates with semigroup theory in extended spaces. In this way obtained a-priori 
estimates allow for a construction of a nonlinear semigroup which evolves finite energy solutions. 



1.7 Statement of Main Results 

Recall, our state space for the analysis to follow is 

Y = H 1 (R^) x L 2 (K+) x H%(Q) x L 2 (Q). 

Theorem 1.1 (Linear). Consider linear problem in fl.6\) with /(it) = 0. Let T > 0. Then, 
for every initial datum (0Qj 4>i] u o, u i) € Y there exists unique generalized solution 



(W),Mt)Mt)Mt))eC([o,T\,Y) 

which defines a Co-semigroup T t : Y — > Y associated with {1. 7| ) (where f = 0). 
For any initial data in 



(1.13) 



Y 1 



v = (</>,0i;u,«i) e Y 



= -[tti + UO x u] ■ l n G H^R 2 ), 
A^ + Ujid^} G La(fi) 



(1.14) 



the corresponding solution is also strong. 



Theorem 1.2 (Nonlinear). Let T > and let f(u) be any nonlinear internal force (Kirchhoff, 
von Karman, or Berger) given by Assumption ] 1. 11 Then, for every initial data (fa, 4>i;uq, u±) € 
Y there exists unique generalized solution (0, <fi t ; u, u t ) to il.6\) possessing property (|1.13p . This 
solution is also weak and generates a nonlinear continuous semigroup St : Y — > Y associated 
with {TTp. 
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If ((j>o, 4>i;uo,ui) € Y\, where Y\ C Y is given by (|1.14|) . then the corresponding solution is 
also strong. 

Remark 1.3. In comparing the results obtained with a subsonic case, there are two major 
differences at the qualitative level: 

1. Regularity of strong solutions obtained in the subsonic case [17) coincides with regularity 
expected for classical solutions. In the supersonic case, there is a loss of differentiability in 
the flow in the tangential x direction, which then propagates to the loss of differentiability 
in the structural variable u. 

2. In the subsonic case one shows that the resulting semigroup is bounded in time, see [161 
Proposition 6.5.7] and also jTTJ 03] • This property could not be shown in this analysis, 
and most likely does not hold. The leack of the energy in energy relation can not be 
compensated for by the nonlinear terms (unlike the subsonic case). 

1.8 Proof Strategy 

In order to orient and guide the reader through various stages of the proof, we briefly outline 
the main steps. 

1. As motivated by the linear theory in the subsonic case, we use the modified energy (as 
given in the previous section) to setup the linear problem abstractly. We decompose the 
linear dynamics into a dissipative piece A (unboxed below) and a perturbation piece P 
(boxed below): 

'{d t + ud x )ct> = ii> 

(d t + Ud x )i> = Atp - [i0 
< d u (j) = -d t u ■ 1q(x) - U 
utt + A 2 u = j[ip] 
u = d v u = 

We then proceed to show that A (corresponding to the unboxed dynamics above) is m- 
dissipative on the state space. While "dissipativity" is natural and built in within the 
structure of the problem, the difficulty encountered is in establishing the maximality prop- 
erty for the generator. The analysis of the the resolvent operator is no longer reducible to 
strong elliptic theory (unlike the classical wave equation) . The "loss of cllipticity" prevents 
us from using the known tools. To handle this we develop a non-trivial approximation 
argument to justify the formal calculus; maximality will then be achieved by constructing 
a suitable bilinear form to which a version of Lax-Milgram argument applies. 

2. To handle the "perturbation" of the dynamics, P (boxed) we cast the problem into an 
abstract boundary control framework. In order to achieve this, a critical ingredient in the 
proof is demonstrating "hidden" boundary regularity for the acceleration potential ip of 
the flow. It will be shown that this component is an element of a negative Sobolcv space 
£2(0, T; iJ -1 / 2 (f2)). The above regularity allows us to show that the term < u x ^[iji\ > 
is well-defined via duality. Consequently, the problem with the "perturbation" of the 
dynamics P, can be recast as an abstract boundary control problem with appropriate 
continuity properties of the control-to-state maps. 



111 



-+ 



x (0,T), 



T d x u ■ lfi(x) 



in W + x (0,T), 

on R {W)} x (°' r )< 



(1.15) 



in n x (0,T), 
in dfl x (0,T). 
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3. Then, to piece the operators together as A + P, we make use of variation of parameters 
with respect to the generation property of A and appropriate dual setting. This yields 
an integral equation on the state space (interpreted via duality) which must be formally 
justified in our abstract framework. This is critically dependent upon point (2.) above. We 
then run a fixed point argument on the appropriate space to achieve a local-time solution 
for the fully linear Cauchy problem representing formally the evolution y t = (A + P)y G Y . 
In order to identify its generator, we apply Ball's theorem [4] which then yields global 
solutions. 

4. Lastly, to move to the nonlinear problem, we follow the standard track of writing the 
nonlinearity as a locally Lipschitz perturbation on the state space Y. In the most difficult 
case of the von Karman nonlinearity the latter is possible due to the established earlier 
"sharp" regularity of Airy's stress function |T6j| pp. 44-45]. This allows us to implement 
local theory with a priori bounds (for T fixed) on the solution. The global a priori bounds 
result from an appropriate compactness-uniqueness argument supported by a maximum 

1 /2 

principle applied to Monge-Ampere equation; the latter implies HmHl^si) < C7| [ [tt, uJHj/ 
for u e H 2 {Q) n H^(Q), see QH Sections 1.4 and 1.5, pp. 38-58] for details. The above 
procedure yields a nonlinear semigroup which, unlike the case of subsonic flow, is not 
necessarily bounded in time (see the case of subsonic dynamics [HI \T7\ [43] ) and the 
resolvent of this semigroup is not compact. 



2 Abstract Setup 



2.1 Operators and Spaces 

Define the operator A = — A + jx with some /i > (H and with the domain 

9{A) ={m£ H 2 {R 3 + ) : d v u = 0}. 

Then ^{A 1 ' 2 ) = H^Rl) (in the sense of topological equivalence). We also introduce the 
standard linear plate operator with clamped boundary conditions: srf = A 2 with the domain 

9{af) = {u e H\n) : u\ dn = d v u\dn = 0} = (H 4 n H 2 )(fl). 

Additionally, &{^ x / 2 ) = (H 2 n H^)(tt). Take our state variable to be 

y = (0, V; u, v) e (^(A 1 / 2 ) x L 2 (K+)) X (^(,t/ 1/2 ) x L 2 (fi)) = Y. 

We work with ijj as an independent state variable, i.e., we are not explicitly taking %p — 4>t + U<j) x 
here. 

To build our abstract model, let us define the operator A : 2>(K) C Y — s- Y by 







( -ud x <f> + i/j \ 


i> 




-Ud x if>- A((f> + Nv) 


u 




V 


\v) 




^ -stu + N*Atp J 



(2.1) 



5 We include the term fil in the operator A to avoid a zero point in the spectrum. After we produce our 
semigroup analysis, we will negate this term in the abstract formulation of the problem in order to maintain the 
equivalence of the abstract problem and the problem as given in 111.71 1. 
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where the Neumann map N is defined as follows: 

Nf = g <^=> (-A + n)g = in R 3 + and d v g = f for z = 0. 

Properties of this map on bounded domains and are well-known (see, e.g., j3Ql P-195] and 
QH Chapter 6]), including the facts that N : H S (R 2 ) i-> H S+1 / 2 (R\) and 

N*Af = 7 [f] for f€9(A), 

and via density, this formula holds for all / <E 9 (A 1 ' 2 ) as well. Additionally, when we write Nv 
for v : CI — > R, we implicitly mean Nv ex t, where v ex t is the extension^ by outside CI. 
The domain of 9(A) is given by 



9(A) 











u 


£ Y 









-ud x (j) + ^£H x {m? + ), 

-Ud x ip - A((t> + Nv) £ L 2 (R%), } (2.2) 

v £ 9(^ 1 / 2 ) = H$(Cl), - srfu + N*A4> £ L 2 (Cl) 



We can further characterize the domain: 

Since on 9(A) we have that ip = Ud x (j)+h for some h £ H l {R\), then 7^] £ i?" 1/2 (R 2 ) (we 

identify R 2 and dR 3 + ). This implies ~/[ip]\ n £ R- 1 ' 2 ^) = [H^ 2 (tt)}' . Therefore we have that 

siu £ H- 1 ' 2 ^) C [9(.<^ x ^)\ (recall that by interpolation the relation 9(<s/ 1 / 8 ) C h] /2 (CI) 
holds). Thus 

u £ 9(^ 7 / 8 ) C H 7 ' 2 {n). 
Moreover, for smooth functions ip £ L 2 (R 3 I _) we have that 

(Ud x ip + = (Ud x ip + A((j) + Nv),^) L2[R3+) - <u,7[^] >l 2 (r2) ■ 

Thus on the account that (<f>,ip;u,v) £ 9(A), so that Ud x ijj ~ A(<f> + Nv) £ L 2 (R\_) and 
v £ Hq(Q) we have that 

\{Ud x ^ + A0,$) R » | < C||^|| R 3 + ||7Mlk-3( K =) 

for any $ £ L 2 (R 3 + ) with 7^] G H- 2 (R 2 ). 

Writing Acf) = (—Ud x ip) + l 2 for some l 2 £ L 2 {R 3 + ) with the boundary conditions d v <j> = v, 
where v £ Hq(CI), wc easily conclude from standard elliptic theory that 

(j) = -UA- l d x ip + h 2 for some h 2 £ H 2 (R\). (2.3) 

Substituting this relation into the first condition characterizing the domain in (|2.2[) we obtain 

U 2 d x A' 1 d x ^ +ip = h!£ H 1 ^) 

which implies 

U 2 d x 2 A- 1 d x ^ + d x ^£L 2 {R\) 

Introducing the variable p = A~ l d x ip one can see that p satisfies wave equation in the supersonic 
case 

(U 2 - \)3 2 p + (-A y . z + fx)p £ L 2 (R 3 + ), (2.4) 

where d v p = on the boundary z = distributionally. 

The observations above lead to the following description of the domain 



5 We must utilize this zero extension owing to the structure of the boundary condition for dv<§> in 111. 151 1. 
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Lemma 2.1. The domain of A, £F(A) C Y, is characterized by: y € S>{A) implies 
. !/=(W,u,i))ey, 7[V>] G ff- 1/2 (^), 

• -Ud x (j) + i\) e H^Rl), 

• we ^K 1 / 2 ) = # 2 (fi), u g ^K 7 / 8 ), 

• |(-t/a,V-^,^) R -;| < oo, V^eL 2 (R^) unft 7 Mer 2 (K 2 ), 

• U 2 d x 2 A- 1 d x 'ip + d x tp e £ 2 (K+) or (OOP Since 6y (JO]) = -Up + h 2 for some 
/i2 S i? 2 (R+), equation can be also written explicitly in terms of <p as 

(U 2 - l)3 x 2 + (-A V!Z + n)(j> € L 2 (R 3 + ) 

where d u (f> = v ex t on the boundary 2 = 0. 

2.2 Cauchy Problem and Unbounded Perturbation in Extended Space 

With this setup, we will be in a position to show that A is m-dissipativc. A peculiar feature 
introduced by the presence of the supersonic parameter is the loss of uniform ellipticity in 
the static version of the perturbed wave operator and the loss of compactness in the resolvent 
operator. The domain of A does not posses sufficient regularity. To cope with this difficulty, 
suitable approximation of the domain will be introduced. In view of this, the proof of the 
maximality property is involved here. The obtained result will give that the Cauchy problem 

yt = Ay, 2/(0) = y a e Y (2.5) 

is well-posed on Y. We will then consider the (semigroup) perturbation 



P 







< \ 


V> 




-UANd x u 


u 







W 




V. J 



The issue here is the unbounded "perturbation" , which does not reside in the state space Y . 
Note that £%{AN} <f_ L^f^A, and only the trivial element is in the domain of AN when the 
latter considered with the values in Y. This fact forces us to construct a perturbation theory 
which operates in extended (dual) spaces. This step will rely critically on "hidden" boundary 
regularity of the acceleration potential ip - established in the next section. As a consequence, 
we show that the resulting Cauchy problem y t = (A + F)y, y(0) = yo <G Y yields well-posedness 
for the full system in (| 1 . 1 5[) . Application of Ball's theorem fj] allows us to conclude that A + P, 
with an appropriately defined domain, is a generator of a strongly continuous semigroup on Y. 

3 Proof of Main Result 

3.1 Hidden Regularity of 7^] 

We consider the following initial boundary value problem: 

(d t + Ud x ) 2 4> = A<f> in R3 , 

0(0) =0 O ; <M0)=</u, (3.1) 
d v (f> = h(x,t) onR 2 {(xy)} . 
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We assume 

h(x,t) G 4 oc (0,T;i 2 (R 2 )). (3.2) 

We note that initially, in our studies of the partial dynamics associated to the dissipative part of 
the dynamics (the system in (|1.15j) with the boxed term removed), we will take h(x, t) — —dtu, 
where u is a plate component of a generalized solution to the unboxed part of (] 1 . 1 5|) . Later, in 
considering the perturbation of the dissipative dynamics, we will take 

h(x, t) = - [d t u + Ud x u] exf , 

where u is some other function which belongs the same smoothness class as u. Therefore the 
requirement in (|3.2j) is reasonable. 

Let (j> be the energy type solution of (|3.1[) . i.e. 

(0, 4> t ) G L 3 (0, T; H 1 (R 3 ^) x L 2 (R 3 + )), VT > 0. (3.3) 

These solutions exists, at least for sufficiently smooth h (see [33J and |38|). 

Our goal is to estimate the trace of the acceleration potential 4>t + U<f> x on z = 0. The a priori 
regularity of <j)(t) G H 1 (M.^.) implies via trace theory <^ x (*)|z=o G i?~ 1//2 (M 2 ). However, the a 
priori regularity of <f>t does not allow to infer, via trace theory, any notion of a trace. Fortunately 
we will be able to show that this trace does exist as a distribution and can be measured in a 
negative Sobolev space. The corresponding result reads: 

Lemma 3.1 (Flow Trace Regularity). Let (|3.2p be in force. If <j)(x,t) satisfies (|3.ip and 
(|3~3)l . then 

dtM, d x j[<p] G L 2 (0, T; ff~ 1//2 (R 2 )) VT>0. 
Moreover, with ip = <f>-t + U(f> x we have 

^ h\m)T H -v^w)<tt < °t \E fl {0) + £ \\d„ct>(t)\\ 2 dtj (3.4) 

The above result is critical for the arguments in later portion of this treatment. Specifically, 
the above result holds for any flow solver; we will be applying this result in the case where 
d v 4> = —v G C(0,T;i2(fi)) coming from a semigroup solution generated by A. 

Proof. One can see that the function rj(x, t) = 0(x + Ute\,t) = <p(x + Ut,y, z, t) possesses the 
same properties as in (|3.3I) and solves the problem 

'd?r/ = Ari in R\, 

r?(0) = <£o; r)t(0) = <t>i + Ud x <j> , 



d v 4> = h*(x,t) on 



p2 



where h*(x,t) = h{x + Ut ei ,t) which is also belong to L 2 (0,T; L 2 (R 2 )), VT > 0. 

Our goal is to estimate the time derivative r\ t on z = 0. Clearly, the a priori regularity of r\ t 
does not allow to infer any notion of trace. However, we will be able to show that 

T \\l[m]{t)\\ 2 H -u^) dt < Ct + \\Vi\\ 2 + £ \\h*(t)\\l 2m dt\ ■ (3.5) 

Since ?/t(x, t) = ip(x + Utei, t), we can make inverse change of variable and obtain the statement 
of Lemma 13.11 
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In order to prove (|3.5j) we apply a principle of superposition with respect to the initial and 
boundary data. 

Step 1: Consider h* = 0. Here we apply Theorem 3 in [33J with k = 0. This yields 

\\v(t)\\Hi { Rl) + \\Vt{t)\\ 2 + ll7Mll| 2(0 ,T;H-i/2(K=)) + I It[»|]| \l 2 (0,T;H^m) 



< C T 



\v(0)\\hhr%) + \M0)\\* ■ (3-6) 



Step 2: Consider zero initial data and arbitrary h* <S L2((0, T) x IR 2 ). We claim that the 
following estimate holds: 

ll7Mt|lL(o,'zw-V2(R2) + ll7MllL(o,T;ffV2(R2)) < Ct[||/i*||I 2 (e)]- (3-7) 

The proof of the estimate in (|3.7[) depends on the fact that the problem is defined on a half-space. 
Since </>0; </>i = we then take the Fourier-Laplace transform: 

t^>T = £ + ia, (x,y) ->i(j, = i(fj, 1 ,H2), 

with £ fixed and sufficiently large. 

Now, let rj = rf(z, (i, a) be the Fourier-Laplace transform of r\ in x, y and t, i.e. 



1 



(2^) 2 

This yields the equation 



K 2 



+00 



dxdy / dte~ Tt ■ e ~^^+v^) . v (x, y, z, t), Rer > 0. 



%* = (M 2 + t 2 -o*+2i£a)fi, 
with transformed boundary condition, 

%(z = 0) = -h*(p,r), T = £ + ia. 
Solving the ODE in z and choosing the solution which decays as z — > +00, we have 

1 



with 



V( z : M) r ) = -7=h*(n,T)exp{-zy/s), 



where the square root y/s is chosen such that Rey/r 2 + |^| 2 > when Re r > 0. On the 
boundary z = we have 77(2: = 0, fj,, r) = -4=/i* (/i,r). Taking the time derivative amounts to 
premultiplying by t = £ + icr, and with a slight abuse of notation, we have 



%{z = 0,/i,r) = —=}i*(ij,,t), t = £ + i<t. 

\/S 



Denoting the multiplier - 



m{a, /z), we can infer the trace regularity of 7/ from m(a, /z): 



(N 2 ) 1/4 



((H 2 - ^ 2 ) 2 + + 2^ 2 + 2|/i| 2 e 2 ) 



1/4 
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Lemma 3.2. We have the following estimate 



\m(t,a,ii)\<2 



1 e 



VK;v)et 4 , e^o. (3.8) 



Proof. First, we can take each of the arguments £, a and \/j,\ to be positive. Next we consider 
the partition of the first quadrant of the (|/z|, a)-plane as follows: 

'(a) \n\<<r/y/2, 

(b) \fi\ > V2a, 

(c) o-/V2 < \fx\ < V2a. 

In cases (a) and (b) above, we can write 



\ m &°>l*)\ ^ ~ " TT7I ^ 7 77/4 ^ 2 - 

(|M a -* a | +e 4 ) (^ 4 /4 + e 4 ) 



In case (c) we have 



i « ^ V^i^lf ^ i V^+e 



(C 4 + 2£ 2 <7 2 + 2H 2 £ 2 ) 1/4 l?l 1/2 (£ 2 + 2a 2 + 2| M | 2 ) 1/4 

]_ a/ct 2 + j 2 1 , 2 2 ^i/4 _J_/ 9 |„|2 , ta\V4 



This implies the estimate in (|3.8|) . □ 
By the inverse Fourier-Laplace transform we have that 

e-^T] t {x,y,z = 0,t) = ^- ( dpi [ dae lat ■ e^ xf " 1+y ^ ■ m(^a, n)'f?(^^ + ia). 

27T ./D2 Lno 



Thus by the Parseval equality 
1 / ,+oc 

^Jo h ~ %{z ~ ""-'"x-^w)- ~ j m2 (i + imp; 

*?t) < (1 + r 2 ) 1/2 / rfM /" rfcr|/?( M , £ + ia)\ 2 



and hence 



Since 



h*(fj,,^ + ia) = —-^ / rfxdy / (ft e~ ttTt • e -«(*/*i+wa) . e -^h*(x,y,t), £ > 0, 



(C;^)<2^(l+r 2 ) 1/2 / rfardy / ttt e-^*|A*(a: f !/,*)| 



we obtain that 



This implies the estimate for ||7[ 7 ?]t|||,( t-h- 1 / 2 ^ 2 ) m P-^P - To obtain the corresponding bound 
for ||7[?7]|| 2 2 ( T-ffi/2(i{2)) we use a similar argument. 

The relations in (|3.6[) and (|3.7[) yield the conclusion of Lemma 13.11 □ 
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Remark 3.1. In Lemma 13.11 we could also take an arbitrary smooth domain O instead of K+. 
Indeed let Q — O x (0,T) and £ = dO x (0, T). Assuming a priori H 1 (Q) regularity of the 
solution, then one can show that j[r)t] £ £2(0, T\ H~ 1 / 2 {dO)) for all smooth domains. The a 
priori H l (Q) regularity is automatically satisfied when the Neumann datum h* is zero and the 
initial data are of finite energy (H 1 x L 2 ). In the case when h* is an arbitrary element of L 2 (T I ), 
the corresponding estimate takes the form 

ll7Mli 2 ( ,T;i?-i/=(]R=)) + ll7Mlli 2 ( ,T;ffi/2(R2)) < Ct[||7i*||! 2( e) + \\v\\m(Q)]- (3-9) 

The proof of this estimate can be obtained via microlocal analysis by adopting the argument 
given in [29l [33] . In our case when O = we have the estimate in (|3.7[) which does not contain 
the term ||??||#i(q) and hence it can be extended to less regular solutions. In the case of general 

domains £2(2) Neumann boundary data produce in wave dynamics only H 2 ^{Q) solutions with 
less regular (than (|3.9[l ) boundary traces 7 [77] £ i? 1 / 3 (E) and 7(774] £ i?~ 2 / 3 (£) (see [32] and 
also |29j ) . The above result is optimal and can not be improved, unless special geometries for 
n are considered. We also note that the above mentioned result improves upon [33] where the 
interior regularity with L 2 (T,) Neumann data is only 77 <G i/ 1 / 2 (Q), rather than H 2 / 3 (Q). Thus, 
for general domains we observe an additional loss, with respect to ()3.7|) . of smoothness for the 
boundary traces 7(77] (1/6 = 2/3— 1/2 of the derivative). 

3.2 Linear Generation of Unperturbed Dynamics 

Our main result in this section is the following assertion: 

Proposition 3.3. The operator A given by (|2.1[) and (|2.2[) is maximal, dissipative and skew- 
adjoint (i.e. A* = —A). Thus by the Lumer- Phillips theorem (see '351) A generates a strongly 
continuous isometry group e At in Y . 

Our calculations for the proof requires some approximation of the domain S>{K) as a pre- 
liminary step. 

3.2.1 Domain Approximation 

We would like to build a family of approximants which allows us to justify the formal calculus 
occuring in the subsequent dissipativity and maximality considerations. Below we concentrate 
on the more difficult supersonic case U > 1. 

Let us take arbitrary (cf),ip;u,v) in S?(A). Then 

-U<j> x +i)&H l {W\), -Uda.il>- A{<f> + Nv) eL 2 (R 3 + ). 

Since we also have that <f> e i? 1 (IR 3 f .), there exist h £ ff 1 (R 3 f ) and g £ (depending on 

the element in the domain) such that 

-U(t> x + tp-rct) = h£H 1 {m? + ), -Ud x ip-A((l) + Nv)=g£L 2 (R 3 h ) (3.10) 

for every rgl. From the first relation: 

d x tp = Ud x 2 <j) + rd x (j> + h x . 

Substituting into the second yields 

~U 2 d x 2 (f> - Urd x <t> - Uh x - A{<j> + Nv) = g. 



1G 



This is equivalent to 

(U 2 - l)d x 2 (ft + Ur<ft x - (A y , z + n)<j> = -g - Uh x = f G L 2 {M. 3 + ) (3.11) 

with the boundary conditions d v (ft — V ex t, where v G Hq(£1). By the trace theorem there exists 
77 G i? 7 / 2 (R^_) such that d v r\ = v ext . Therefore we can represent 

(ft = (ft. t + T), 

where </>* G i/ 1 (K.^_) is a (variational) solution to the problem 

(U 2 - l)0 x 2 ^ + Ur<j>* x - {A y , z + fi)4>, = f, G L 2 (H5-), (3.12) 
d v cf>* = 0, 

and 

f*=f-(U 2 - l)d x 2 r] - Urn x + (A y , z + /i)r/. 

Due to zero Neumann conditions the problem in (|3.12|) can be extended to the same equation in 
R 3 . Therefore we can apply Fourier transform in all variables. This gives us (x u), (y, z) <-> 
k G M. 2 ) : 

( C{juj 2 + irUto + |fc| 2 + (i)<ft* = /* G L2, 
where = (U 2 — 1). Sinccfl 

I - cjju 2 + iUru + \k\ 2 + A i| 2 = r 2 U 2 uo 2 + (\k\ 2 +fi- c v uj 2 ) 2 > c {\k\ 2 + uj 2 + 1) 

with some cq > for appropriate r = r(cu, £t) > 0, the formula above leads to the solutions <ft* in 
i? 1 (]R 3 f ) for every /* G L 2 (R+)- Moreover if /* G iP(R 3 ) is such that its even (in z) extension 
on M 3 has the same smoothness, then 0„ G i2~ 1+s (R+) for every s > 0; hence for s < 3/2 we can 
define a continuous affine map 

t : H"(R%) -> H S+1 (R% ), where r(/) = 0, 

where ft solves (|3.11[) with the Neumann boundary condition d v <j> — v ext G H 2 (M. 2 ). 

Thus, in order to find an approximate domain T> n it suffices to solve (|3.11[) with the right 
hand side in _ff 1 (R 3 f ). Hence, we arc looking for <ft n G i? 2 (R 3 f ) such that 

(U 2 - \)d 2 $ n - A y . z <ft n + ficj) n + Urft n x = -g n - Uh n x G H\R 3 + ), (3.13) 

and d v <ft n = v ex t for all n, and where h n — > h in i/ 1 (R 3 f ) and g n — > g in ^(R 3 ^). Here h and g 
are given by (|3.10[) . 

Solving equation (|3 . 1 3[) with right hand side — g n — Uh x G i? 1 (R 3 f ) gives solution 

(t> n e H 2 {R 3 + ), <p x G iJ 1 (R 3 f ). 

We then define V" = U<j% + r<ft n + h n G i? 1 (R 3 f ). Then 

Ui>l + A{cft n + Nv) = U 2 (ft n xx + Uh n x + A((ft n + Nv) + rUtfe = -g n -> g = Uip x + (-A + n)<j> 

as desired. Additionally, since r is affine and bounded, and each <ft n has the same boundary 
condition, we can conclude that 

110" - cft\\ m(K) = \\r(g n + UK) - r(g + Uh x )\\ m(R3+) < \\ft n - <ft + U(h" x - MIl,^) ^ «■ 

We may proceed similarly on ip n = U<ft x + r<ft n + h n . Thus, we have obtained 

7 In the subsonic case (cjy < 0) the equation in 1 13.1211 is elliptic and we can obtain better estimate. 
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Lemma 3.4 (Domain Approximation). For any y = (<f>,ip;u,v) G 3(A) there exist approx- 
i\),i> n G i/ 1 (R^_) 



imants cf> n G H 2 (R\),ip n G H 1 (Ri) swcft tfmt y n = ((j) n ,ip n ;u,v) G 9(A) and y n -> y in Y . 



Moreover 



As a consequence 



Uip% + A{cj) n + Nv) -> Uip x + A{<f> + Nv), in L 2 (R+), 
V -U<f%->i/>- U<j> x , in H 1 ^ 3 ,). 



(Ay n ,y n ) Y (Ay,y) Y for all ye 3(A). 



3.2.2 Dissipativity 

The above approximation Lemma allows us to perform calculations on smooth functions. 

Let y n G 31(A) be the sequence of approximants as in Lemma [3.41 First, we perform the 
dissipativity calculation on these approximants (which allows us to move A 1 ' 2 freely on flow 
terms <j> n and ip n ): 



( 


-Ud x (f> n + ip n 




' (j) n 


\ 




-Ud x ip n ~ A(4> n + Nv) 










V 




u 




V 


-£/u + N*Ai/j n 




V 


I 



(Ay n ,y n ) Y = 

U Lb 
_ rJn, _L AT* AJ,n 

@(A 1 / 2 )xL 2 ((l)x3>(stf) 1 / 2 xL 2 (n) 

= (A 1 / 2 (V' n - Ud x (j) n ),A 1 / 2 (j} n ) - (Ud x ^ n + A((j) n + Nv),ip n ) 

+ < £/ 1/2 v, £/ 1/2 u >-<srfu- N*A^ n ,v > 
= - U(A 1 f 2 d x (t> n ,A 1 / 2 (t> n ) - (Ud x ip n + ANv,iP n )+ < N*Aip n ,v> . 

One can see that 

(A 1 / 2 d x ab n ,A l l 2 ah n )= [ Vd x r ■ Mf 1 = I [ d x \V^\ 2 = 

Similarly (Ud x tp n ,^ n ) = 0. Therefore 

(Ay™, y n ) Y = -< ANv, ip n > + < N*Aip n ,v > = 0, 
Furthermore, by the convergence result in Lemma f3.4[ we have that for all y G 3(A) 

(Ay,y) Y = lim (Ay n ,y n ) Y = 0. 

n— >oo 

This gives that both operators A and —A are dissipative. 
3.2.3 Maximality 

In this section we prove the maximality of the operators A and —A. For this it is sufficient to 
show M{X -A)=Y for every A £ t, A ^ 0, i.e. for a given F = ((/>', ip'; u', v'), find a V G 9(A) 
such that (A — A)V = F. Writing this as a system, we have 



'X^ + Ud^-tp = 0' G 9(A 1 / 2 ) 1 

\i) + Ud x ip + A((f> + Nv) =i(>' €L 2 (R 3 + ), 

\u-v = u' G 3(£? 1 ' 2 ), 

Xv + srfu - N*Aip =»'ei 2 (f!), 



(3.14) 
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(recalling that Nv is taken to mean Nv ex t where v ex t is the extension by zero outside of f2). 
In the space Y we rewrite (|3.14j) in the form 

a(V,V) = (F,V) Y , (3.15) 

where for V = {<fi, ip; u, v) and V = (4>, ip; u, v) we denote 

a(V, V) =(X(f> + Ud x( f> - ^ 

+ (Xip + Ud x i> + A{<j) + Nv), i/i) R 3 

+ (Xu - v, jrfu)n + (Aw + srfu - N*Aip, v)q. 

Let {r/k} x {efe} be a sufficiently smooth basis in ^(A 1 / 2 ) x &(s^ x / 2 ). We define an iV- 
approximatc solution to p.ISp as an element 

V N e Y N = Span {(r) k , rjf,e m , e„) : 1 < k, I, m, n < N} 

satisfying the relation 

a(V N ,V) = (F,V) Y , WeY N . (3.16) 

This can be written as a linear 4N x AN algebraic equation. Calculations on (smooth) elements 
V from Ijv gives 

a(Y,V) - A{||^|| 2 + i^ll 2 + \W 1/2 u\\ + imi 2 } 

This implies that for every A =^ the matrix which corresponds to (|3.16|) is non-degenerate, and 
therefore there exists a unique approximate solution Vjv = {4> N , ^ JV ; u N ,v N ) . Moreover we have 
that 

a{V N ,V N ) = (F,V N ) Y 

which implies the a priori estimate 

P 172 ^ + |h^|| 2 a + \W 1/2 U N \\l + \\V N \\1 < ^\\Ffy 

Thus {V)v} is weakly compact in Y. This allows us to make limit transition in (|3.16p to obtain 
the equality 

a{V,V) = {F,V) Y , VV£Y M , VM, 

for some V £ Y. Thus p.!4p is satisfied in the sense distributions. This proves maximality of 
both operators A and —A. 

Since both operators A and —A are maximal and dissipative, we can apply [T2J Corollary 
2.4.11. p. 25] and conclude that the operator A is skew-adjoint with respect to Y. This completes 
the proof of Proposition 13.31 

The fact that A is skew-adjoint simplifies calculations later in the treatment. In what follows, 
we use 3>(K) and fi?(A*) interchangeably. 

Remark 3.2. To conclude this section, we mention that for yo £ Y the Co-group e At generates 
a mild solution y(t) = e At yo to the PDE problem given in (jl.151) without the boxed term (see 
also (|3.17[) below with k = 0). We also note that adding a linearly bounded perturbation to the 
dynamics will not affect the generation of a Cg group. In particular, the addition of internal 
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damping for the plate of the form kut, k > on the LHS of does not affect generation and 
we can construct Co-group Tk(t) which corresponds to the problem 

' {d t + Ud x )(j) = tp 

< d y (j) = -u t ■ ln(x) 
u t t + ku t + A 2 u = j[ip] 
u = d v u = 

Moreover, this damping term ku t does not alter £>(A) or 3>(A*), and hence the analysis to 
follow concerning the full dynamics (A + P, and the addition of nonlinearity) is valid presence 
of interior damping. We plan use this observation in future studies of this model. 

3.3 Variation of Parameters and Perturbed Linear Dynamics 

We begin with some preparations. 

3.3.1 Preliminaries 

We would like to introduce a perturbation to the operator A which will produce the non- 
monotone flow-structure problem above. For this we define an operator P : Y — > &(P) as 



(<t>\ 




( \ 






-UANd x u 


u 


= p#M = 





w 




1 o / 



(3.18) 



Specifically the problem in (|1.15[) has the abstract Cauchy formulation: 

Vt = (A + F)y, y(0) = y , 

where ya € Y will produce semigroup (mild) solutions to the corresponding integral equation, 
and yo £ S)(K) will produce classical solutions. To find solutions to this problem, we will consider 
a fixed point argument, which necessitates interpreting and solving the following inhomogeneous 
problem, and then producing the corresponding estimate on the solution: 

y t = Ay + P # u, t > 0, y(0) = yo, (3.19) 

for a given u. To do so, we must understand how P acts on Y (and thus P# on Hy(£l)). 

To motivate the following discussion, consider for y £ Y and z = (4>, ip; u, v) the formal 
calculus (with Y as the pivot space) 

(Fy,z) Y =(P # M,z) y = -U{ANd x u,4>) = -U < d x u,i$\ > . (3.20) 

Hence, interpreting the operator P (via duality) is contingent upon the ability to make sense of 
7[-0], which can be done if 7^] £ i7 _1 / 2 (il). In what follows, we show that the trace estimate 
on tp for mild solutions of (|3.19[) allows us to justify the formal energy method (multiplication 
of (|3 . 1 9(1 by the solution y) in order to perform a fixed point argument. 

To truly get to the heart of this matter, we must interpret the following variation of param- 
eters statement for u £ CQR+; Hq (0)) (which will ultimately be the solution to (|3.19)1 ): 

y(t) = e M y Q + [ e A ( 4 ^)p # [u(s)]ds. (3.21) 
Jo 



in R 3 + x (0,T), 
in R 3 + x (0,T), 

onM ky)} x (°' r )' 
in Q. x (0,T), 

in dn x (0,T). 



(3.17) 
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To do so, we make use of the work in [30] and write (with some AeR, A 7^ 0): 

y(t) = e At y + (A - A) f - A)- 1 P # [n(s)]ds, (3.22) 

Jo 

initially interpreting this solution as an element of [S#(A*))' = [£F(A)]', i.e., by considering the 
solution y(t) in (|3.22[) above acting on an element of f^(A*). 

3.3.2 Abstract Semigroup Convolution 

At this point we cast the discussion of the perturbation P (acting outside of Y) into the context 
of abstract boundary control. By doing this we simplify and distill our exposition, and moreover, 
provide a context for further boundary control considerations. For this discussion, we select and 
cite some results from [301 pp. 645-653] which will be used in this section. 
Let X and U be reflexive Banach spaces. We assume that 

(CI) A is a linear operator which generates a strongly continuous semigroup e At on X . 

(C2) B is a linear continuous operator from U to )]' (duality with respect to the pivot 

space X), or equivalently, (A - A)~ 1 B g 3?(U,X) for all A 6 p(A). 

For fixed < T < 00 and u S Li(0, T; U) we define the convolution operator 

(Lu)(t) = [ e A(t - s) Bu{s) ds, < t < T, 
Jo 

corresponding to the mild solution 

x(t) = e At x + (Lu)(t), 0<t<T, 
of the abstract inhomogeneous equation 

x t = Ax + Bu e [®{A*)]\ x{Q) = x , 
with the input function Bu(t). 

Theorem 3.5 (Inhomogeneous Abstract Equations). Let X and U be reflexive Banach 
spaces and the conditions in (CI) and (C2) be in force. Then 

1. The semigroup e At can be extended to the space [S>{A*)\' . 

2. L is continuous from i p (0,T; X) to C(0,T;[@(A*)]') for every p £ [l,oo]. 

3. Ifue C l (Q,T;X), then Lu g C(0,T;X). 
4- The condition 

(C3) There exists a constant Ct > such that 

/ \\B*e A V|£.di <C T \\x*\\ 2 x ,, Vi* e @(A*) c X*, 
Jo 

is equivalent to the regularity property 

L: £2(0, T; U) — > C(0, T; X) is continuous, 
i.e., there exists a constant kx > such that 

\\Lu\\ C (0,T;X) < kT\\u\\ L2 (o t T;U)- (3.23) 
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5. Lastly, assume additionally that A generates a strongly continuous group e At (e.g., if A is 
skew- adjoint) and suppose that L : L2(0,T;U) — > £2(0, T; X) is continuous. Then (C3) is 
satisfied and thus we have the estimate in (|3.23p in this case. 

3.3.3 Application of the Abstract Scheme 

We now introduce the auxiliary space which will be needed in the proof of the next lemma: 



e Y : -Ud x cl> + i> G H L (R% 



endowed with the norm 

\\y\\z = \\v\\y + II - Ud x 4> + i/j\\ h i {M 3 +) . 

One can see that Z is dense in Y. We also note that by Lemma 12.11 0(A) C Z and thus 
Z' C [0(A)]' with continuous embedding. 

Lemma 3.6. The operator P# given by (|3.18[) is a bounded linear mapping from Hq(Q) into 
Z' . Moreover, the following estimates are in force: 

l|P#MIU' <Cu\\u\\m(n), V u eff 2 (fi), (3.24) 

and also (with A € M, A 7^ 0): 

||(A-A)- 1 P # [ u ]|| y < C UtX \\u\\ H 2 (n) , VueH%(n). (3.25) 

In the latter case we understand (A — A) -1 : [0(A)]' n- Y as the inverse to the operator A — A 
which is extended to a mapping from Y to [0(A)]'. We also have that (pT24|) and (|3~25|) imply 
that P maps Y into Z' and 

\\Vv\\z><C\\y\\ Y and ||(A - A) _1 Py||y < C||y||y, Vy e Y. 

Proof. For u <E H^fl) and y = (0, ip;u,v) <G Z, from p.20p we have 

|(P # M,y)y| = U\(ANd x u,^)\ = U\< d x u,^\ > \. 

Since 7^] = 7[— Ud x <f> + i/j] + Ud x j [<fi], we have from the trace theorem that 

Il7[^]llff- 1 /2 (R 2 ) < C \\(-Ud x ~4> + ip)\\m(s.\) + \U>\\hi(rI) <C\\y\\ z . 



Therefore 



Thus 



*#[u],y) Y \ < C([/)P. T u|| H i/2 (0) ||7[^]|| ff -i/2 (R 2 ) < C(U)\\u\\ H 2 {n) \\y\\ z . 



*#[u]\\z> = sup{|(P # [u] I y)y| :yeZ, \\y\\ z = l} < C(U)\\u\\ H 2 (Q) , 



which implies (|3.24p . The relation in (|3.25p follows from p.24[) and the boundedness of the 
operator (A - A)" 1 : [0(A)]' H> Y. □ 

We may now consider mild solutions to the problem given in (|3.19[) . Applying general results 
on Co-semigroups (sec [35 ) we arrive at the following assertion. 



22 



Proposition 3.7. Let u e C 1 ([0, T];H§(Q)) andy a eY. Then y{t) given by (|3.21j) belongs to 
C([0,T];Y) and is a strong solution to (|3.19|) in [S>{A)]' , i.e. in addition we have that 

and (j3~T9)) holds in [S>(A)]' for each t E (0,T). 

Proof. The result follows from the integration by parts formula 

(Lu)(t)= f e A{t - s) V#[u(s)]ds = (A -A) / e Ait - s \\ - A)- X P # [u(s)]ds 



A / e A(t " s) (A- A)- 1 P # [u(s)]ds- / e A( *" s) (A- A)- 1 P # [-^-?l(.s)]ds 



- e^*(A - A)- 1 P # [u(0)] + (A - A)- 1 P # [u(<)] 
G C(0,T;X), 

and by applying [35J Corollary 2.5, p. 107] to show that y is a strong solution in [£^(A)]'. □ 
Proposition 13.71 implies that y(t) satisfies the variational relation 

d t (y(t),h) Y = -(y(t),Ah) Y + (V#[u{t)],h) Y , Vfte#(A). (3.26) 



In our context in application of Theorem 13.51 we have that X = Y, where Y is Hilbert space 
given by (fI~T2"j) . U = H$(to) and B = P # is defined by ([3~TS]) as an operator from U = H$(Q) 
into Z' (see Lemma I3~6")) . One can see from (|3.20p that the adjoint operator P^ : Z i-> H~ 2 (Vl) 
is given 

P^z = C/[9 :c AfM^]| n = [/a x 7[^| n , for z = 



t 
u 



e 



Condition (C3) in Theorem 13. 5f 4) is then paraphrased by writing 

y i-> P* # e A * {T - t} y = P* # e M y T : continuous Y -> L 2 (0, T; H~ 2 (n)) , 

where ?/t = e A T y = e~ AT y (we use here the fact that A is a skew-adjoint generator). 

Let us denote by e M yr = wr(t) = ((f)(t),ip(t);u(t),v(t)) the solution of the linear problem 
in (|2.5|) with initial data yx- Then from the trace estimate in Q3.4p we have that 



%e M yT\\Uo,T:H-* m = U / \\d x j[m]\\H- H n) 



<C I Hm\\\ 2 H -^ m dt < C t U f i(0) + jf* \\v(t)\\l 2{ n } J dt 

<Ct\\vt\\y = C T ||e- AT y||y = C r ||y||y. 

In the last equality we also use that e At is a Co-group of isometries. 

Now we are fully in a position to use Theorem 13.51 which leads to the following assertion. 



Theorem 3.8 (L Regularity). Let T > be fixed, y Q eYandue C([0, T];Hq(CI)). Ther 
the mild solution 



y(t) = e At y + L[u]{t) = e At y + I ' e A ^-^P # [u(s)}d. 
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to problem (|3.19j) in [S>(A)]' belongs to the class C([0,T];Y) and enjoys the estimate 

max l|y(r)|| Y <||y || y + MNk (0 «n)), Vte [0,T]. (3.27) 

Remark 3.3. The discussion beginning at Theorem 13.51 and ending with the estimate in (|3.27|) 
demonstrates how the perturbation P acting outside of Y is regularized when incorporated into 
the operator L; namely, the variation of parameters operator L is a priori only continuous 
from L 2 {Q,T-U) to C(0,T;[9(A*)]'). However, we have shown that the additional "hidden" 
regularity of the trace of if) for solutions to (|1.7|) allows us to bootstrap L to be continuous from 
£2(0, T;W) to C(0,T;Y) (with corresponding estimate) via the abstract Theorem 13.51 This 
result essentially justifies formal energy methods on the equation (|3.19|) in order to set up a 
fixed point argument (which will follow in the next section). 

For completeness, we also include a direct proof of an estimate of the type (|3.27p in the 
Appendix, independent of the abstract boundary control framework presented in Theorem 13. 51 

3.4 Construction of a Generator 

Let X t = C((0,t];Y). Now, take y = (ej>,lj>; u, v) <E X( and yo <E Y, and introduce the map 
J£" :y —> y given by 

y(t) =e M y + L[u](t), 

i.e. y solves 

y t =Ay + P#u, y(0) =yo, 
in the generalized sense, where P# is defined in (|3.18|) . It follows from (|3.27p that for y 1 , y 2 e X t 

WVi ~ &Vi\Wt < k r\\ui ~u 2 \\ L2{0 j. H 2 {n}} 

< k T Vt max ||ui -U2\\m<m < k T Vt\\y 1 -y 2 l|x t - 

r€[0,t] 

Hence there is < t* < T and q < 1 such that 

ll^yi-^lk <q\\Vi-y 2 h t 

for every t S (0, t*]. This implies that on the interval [0,i*] the problem 

y t =Ay + ¥ yi t> 0, y(0) = y , 

has a local in time unique (mild) solution defined now in Y. This above local solution can be 
extended to a global solution in finitely many steps by linearity. Thus there exists a unique 
function y = (<f), ip; u, v) G C(E+; F) such that 

y(t) = e M yo + f e A ^- s ^[y{s)}ds in Y for all t > 0. (3.28) 
Jo 

It also follows from the analysis above that 

\\y(t)\\ Y <C T \\y a \\Y, t€[0,T], VT>0. 

Thus the problem (|3.28|) generates strongly continuous semigroup T(t) in Y. Additionally, due 
to (|3.26[) we have 

(y(t),h) Y = (y ,h) Y + [ [-(y(r),Ah) Y + (¥[y(T)],h) Y ]dT, V^G^(A), t > 0. 
Jo 
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Using the same idea as in subsonic case [17] (which relies on Ball's Theorem [4] and ideas 
presented in pj5]), we can conclude that the generator A of T(t) has the form 

Az = Az + Pz, z £ 3)(k) = {z £ Y : Az + Fz £ Y} 

(we note that the sum Az + ¥z is well-defined as an element in \3i(A)]' for every z £ Y). Hence, 
the semigroup e At yo is a generalized solution for yo £Y (rcsp. a classical solution for yo £ &(A)) 
to (|TT5j) on [0, T] for all T > 0. 

-Ud x ijj- A((/) + N(v + Ud x u)) £ L 2 (R 3 + ) \ (3.29) 

w g @{sz? 1/2 ) = -*/u + N*Ail) g L 2 (Q) J 

Now we can conclude the proof of Theorem ll.ll in the same way as in |17j by considering bounded 
perturbation of generator of the term 

C(y) = (0,^0,0). 

Indeed, the function y(t) is a generalized solution corresponding to the generator A + P with 
the domain defined in (|3.29|) . This proves the first statement in Theorem ll.il As for the second 
statement (regularity), the invariance of the domain @(A + P) under the flow implies that 
solutions originating in Y\ will remain in C([0,T]; Y\). After identifying <j) t = ip — Ud x <j) one 
translates the membership in the domain into membership in Y\. Elliptic regularity applied to 
biharmonic operator yields the precise regularity results defining strong solutions. The proof 
of Theorem II .11 is thus completed. 



+ P) = { y £ Y 



3.5 Nonlinear Semigroup and Completion of 
the Proof of Theorem 11.21 

To prove Theorem 11.21 we also use the same idea as in [17]. As an intermediate step we obtain 
the following theorem. 

Theorem 3.9. Let ,^(y) = (0, Fi (</>); 0, F 2 {u)) where 

Fi : iJ x (P4) -> i 2 (R+) and F 2 : H$(Q) -t L 2 (0) are locally Lipschitz 
in the sense that every R > there exists cr > such that 

||Fi(0) -Fi(0*)|| R 3 <c R ||^-^|| liE 3 and \\F 2 (u) - Fi(u*)\\ n < c R \\u - u*\\ 2>Q 

for all 4>,<j)* £ H^Rl) and u,u* £ H^{Q) such that ||0|| liR s , ||0*|| 1>R 3 , ||u|| 2 ,n, ||«*||a,n < R- 
Then the equation 

y t = (A + V)y + &{y), y(0) =y £Y (3.30) 

has a unique local-in-time generalized solution y(t) (which is also weak). Moreover, for yo £ 
&(A + F), the corresponding solution is strong. 

In both cases, when t ma x(yo) < oo, we have that ||y(£)||y — > oo as t /* t rnax (yo). 

Proof. This is a direct application of Theorem 1.4 [3S] p. 185] and localized version of Theorem 

1.6 [35] p. 189]. □ 

In order to guarantee global solutions, one must have more information on the nature of the 
nonlinear term. The following result provides relevant abstract conditions imposed on nonlinear 
terms which can be verified: 
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Theorem 3.10. We assume that f is locally Lip schitz from H§(Q) into L 2 (fl) and there exists 
C 1 -functional II(it) on Hq(Q,) such that f is a Frechet derivative ofH(u), f{u) — II' (u). More- 
over we assume that Tl(u) is locally bounded on H^lfl), and there exist r\ < 1/2 and C > such 
that 

r?||Au||^ + n(u) + C> 0, VueH%(n). (3.31) 
Then the generalized solution in i3. 30\) is global, i.e. t max = oo. 

Proof. The relation in (|3.31j) implies that the full energy £{t) defined in (jl.8j) admits the esti- 
mates 

£(t) > c + \\4>\\l % + \\u t \\l + \\Au\\ty - ci = c \\y(t)f Y - ci 

for some positive Cj. Therefore using the energy relation in (|1.9[) on the existence time interval 
and the flow trace regularity (see Lemma l3.ip we can conclude that 

co\\y(t)\\l-<c 1 + £(0)+ f \\Au(T)\\nU(T)\\ H - 1/H1!ll) dT 
Jo 

<C(y ) + C T f \\y(r)\\ 2 Y dr- 
Jo 

Thus the solution y(t) cannot blow up at finite time, i.e., t max = oo. □ 

Thus, in order to complete the proof of the first part of Theorem ll.2l one needs to verify that 
the nonlinear forcing f(u) given in Assumption 1 1 . 1 1 comply with the requirements of Theorem 
[331 

3.5.1 Verification of the Hypotheses 

We note the examples of forcing terms described above satisfy conditions of Theorem 13.101 

Step 1: Kirchhoff model. Indeed, in the case of the Kirchhoff model, the embeddings H S (Q) C 
Loo(ri) for s > 3/2 implies that 

\\f(ui) - f(u 2 )\\n < Cr\\ui - U2||jya-«(fi) (3.32) 

for every U{ E Ho{Q) with ||u||j?2(f2) < R- The functional II(m) has the form 

U(u) = [ F(u{x))dx with F(s) = [ f(£)d£ 
Jn Jo 

It follows from (|1.2p that there exist 7 < Ai and C > such that F(s) > -js 2 /2 - C for all 
sei This implies (pPTj) . 

Step 2: Von Karman model. In the case of the von Karman model the arguments are more 
subtle. We rely on sharp regularity of Airy stress function [23] and also Corollary 1.4.5 in [16]. 

||A" 2 [u,-u;]|| H /2,oc(s- 2) < C||u|| 2 ,ri|Mkn 

where A 2 denotes biharmonic operator with zero clamped boundary conditions. The above 
yields 

\\v(u)\\ W 2,*o {n) < C||u||^ n 
which in turn implies that the Airy stress function v (u) defined in (|1.3jl satisfies the inequality 

||[ui,«(ui)] - [u 2 ,v(u 2 )]\\n < C*(j|uij|2, f2 + ||u|| 2 .o||2)ll u i - u 2 \\ 2 ,n (3.33) 
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(see Corollary 1.4.5 in [E|). Thus, f(u) = — [u,v(u) + Fo] is locally Lipschitz on Hq(CI). 
The potential energy II has the form 

Il(u) = ±Jj\v(u)\ 2 -2([u 1 F })u\ dx 

and possesses the property in (|3.31[) . see, e.g., Lemma 1.5.4 [THl Chapter 1]. It is worthwhile 
to note that the property (|3.31[) is related to the validity of the maximum principle for Monge 
Ampere equations. For functions u G H 2 (fl) one has 

diamO 1/2 

sup u < sup u H -=— 1 1 [u, u\ 1 1 / 

n an V 71 " 

Thus for u G H 2 (Sl) H H^{il) we have (see Lemma 1.5.5 in [IS]): 

diamft 1/2 
maxKaOl < -j= mMU^n) 

The above uniqueness property is critical in proving (|3.3ip for any function u E H 2 (0) n Hq (f2) . 

5tep 3: Berger's model. One can also see that the Berger model satisfies (|3.32|l with S = 
and (|3.31[) holds; for details see [T31 Chapter 4] and [TSJ Chapter 7]. 



3.5.2 Completion of the Proof of Theorem 11.21 

For the proof of Theorem 11.21 it suffices to apply Theorem 13.101 along with the estimates stated 
above. These estimates assert that the hypotheses of Theorem 13.101 have been verified for all 
three nonlinear models under consideration. 

Concerning strong solutions we notice first that on the strength of the estimate (|3.33[) the 
domain of (A + P) + & in the von Karman case is the same as the domain + P). The same 
holds for the other two models. The local Lipschitz property of the nonlinear terms, along with 
global bounds on solutions, allows us to claim the invariance of the domain of the nonlinear 
flow. Thus for the initial data in Y\ one has that the solution y = (<p,^> t ;u,ut) € Loo(0,T; Y\) 
(as in the argument of Theorem 13.91 we refer to Theorem 1.6 [351 p. 189]). This implies 

ye C([0,T];Y), <ft € L^O, T; ff 1 ^)), u t £ L^O, T; H 2 (n)), 
<f> tt = -2cf> xt - U 2 xx +A<f> G £00(0,^X3^)), 

u tt = -A 2 u + 7 [0 t + U4> x ] - /(«) = ~A 2 u + U-yfa] + 7 [0t] - /(«) G Loo(0, T; L 2 (fi)), 
the above implies via elliptic theory and Sobolev's embeddings 
A 2 u G Lo^O, T; ff-V2(fi)) u G ^(0, T; ^ 7 / 2 (r!)), 

(C/ 2 - l)<fe* - ^ - 4> m G ioo(0, T; L 2 (ffit))> ^U=o = + ^ e Loo(0, T; i 2 (M 2 )). 

The above relations imply the regularity properties required from strong solutions. 

The regularity postulated for strong solutions is sufficient in order to define variational 
forms describing the solutions. The existence and uniqueness of weak solutions follow by viewing 
generalized solutions as the strong limits of strong solutions. This, along with Lipschitz estimates 
satisfied by nonlinear forces, allows a passage with the limit on strong solutions. This completes 
the proof of Theorem 11.21 

In conclusion we note that the well-posedness results presented in this treatment are a 
necessary first step in studying long-time behavior of solutions. This can be done without the 
addition of damping mechanisms (see [T3] and [THl Remark 12.4.8]) or in the presence of control 
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theoretic damping, e.g. boundary or interior dissipation (see [HI [32]. In either case, the next 
step will be to show the existence of global attracting sets for the plate component of the model, 
and analyze their properties (i.e., compactness, dimensionality, and regularity). 
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5 Appendix 

5.1 Direct Proof of Estimate (13.271) for Fixed Point Statement 

Letue C* 2 ([0,T];ij2(r>)) and let y(t) = (<j>(t),ip(ty,u(t),v(t)) £ C([0,T];Y) be a a mild solution 
to (|3.19[) . This implies that y(t) is a (distributional) solution to problem 



in R% x (0,T), 
in R 3 + x (0,T), 



(d t + ud x )4> = V 
(d t + Ud x )ij) = A4> - ^ 

d v 4> = -(d t u + Ud x w)-ln{yi) onM| (xi/)} x (0,T), 

w tt + A 2 w = 7 [^] inftx(0,T). 

u = d v u = indfix(0,T). 



(5.1) 



It follows from the trace theorem that there exists r\ from the class C 2 ([0, T]; if 2 (R^)) such that 



duV = —U [d x u] ext on 



(0,T). 



Let <j> = <t>- r). Then it follows from (HD that y(t) = ((f>(t), tp(t); u(t),v(t)) € C([0, T]; Y) solves 
(inhomogeneous) problem 



(d t + ud x )<t> = v + h 

(d t + U8 x )i> = A0 - ficj) + h 
d„4> = -d t u ■ ln(x) 
u tt + A 2 u = 7[?/>] 
u = d v u = 



in R3 x (0,T), 
in R3 x (0,T), 
on ^\ {x>y)} x (0,T), 
in n x (0,T), 
in dfl x (0,T). 



(5.2) 



where 

fi = -(d t + Ud x ) V , / 2 = (A- M )r?. 
Problem (|5 .2[) can be written in the form 

y t =Ay + F(t), y(0) = y , 



(5.3) 



where F = (fi, f 2 \ 0, 0) e ^([O, T]; Y). Therefore by Corollary 2.5[35l p.107] for any y g ^(A) 
there exists a strong solution y to (|5.3[) in Y. This solution possesses the properties 



yeC((0,T);Y)nC 1 ((0,T);^(A)') 
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and satisfies the relation 



\\y(t)\\Y = \W)\\y+ f (F(r),y(r))Ur. (5.4) 
Jo 

Now we can return to the original variable y(t) = (</>(t) = 4>(t) + r)(t),i/j(t);u(t),v(t)) and show 
that (|5.4p can be written in the following way 



= ||tf(0)|& - 2U [\u x (r)^mr)}) L2{n) dr, (5.5) 



provided u £ C 2 ([0, T]; fl^fi)) and y £ 2>(k). 

The integral term in (|5.5[) can be estimated as follows: 







— , 

CO 




\\\yo\\v -t 


/'[ 




Jo L 



2 

Lff 2 (o) 



r/r 



| < u x ,j[ip] > | < Co / [\\w\\ H 2 m dT + \\lbP}\\ H -i/2 m \ dr (5.6) 
J o 

One can see that the estimate in (|3.4[) can be written with the constant Ct which is uniform at 
any interval, i.e. in the form 

hbP](r)\\ 2 H -y 2m dr < C T (e^O) + J \\d^{r)\\ 2 dr 

which holds for every t £ [0, T]. Since in our case d v 4> = — (v + Ud x u) ■ ln(x), we have that 

^ 

for every t £ [0, T]. Therefore (JOJ) and (JSU) yield 

\\y(t)\\ 2 Y < C T (\\yoWy + J Q \HT)\\ 2 H2m dT + j* \\y(r)\\ 2 Y dr 

for every t £ [0,T], where y £ Y\ and and u £ C 2 ([0, T]; Hq (CI)). Now we can extend this 
inequality by continuity for all y £ Y and w £ C([0, T];H$(Cl)) to obtain ([X27]) . 
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